EXPERIMENTS WITH MODULI OF QUADRILATERALS II 



ANTTI RASILA AND MATTI VUORINEN 

Abstract. The numerical performance of the AFEM method of K. Samuelsson 
is studied in the computation of moduh of quadrilaterals. 



1. Introduction 

The moduli of quadrilaterals and rings are some of the fundamental tools in 
geometric function theory, see [Ahl] . |AV Vj . [Kiihj . [LVj . The purpose of this paper 
is to report on our experimental work on the numerical computation of the moduli 
of quadrilaterals, based on the algorithms and software of |BSVj and motivated by 
the geometric considerations in |HVVj and |DVj . The methods considered here may 
be classified into two classes: 

(1) methods based on the definition of the modulus and on conformal mapping 
of the quadrilateral onto a canonical rectangle, 

(2) methods based on the solution of the Dirichlet-Neumann problem for the 
Laplace equation. 

With the exception of a few special cases both methods lead to extensive numer- 
ical computation. For both classes of methods there are several options, see |Gaij . 
[Hen], [Pap]. Among other things, historical remarks are given in [Porj . 

We study the case of a polygonal quadrilateral and the way its modulus depends 
on the shape of the quadrilateral. Following the approach of |BSV] our main method 
is the adaptive finite element method AFEM of Klas Samuelsson and it belongs to 
class (2). We compare this method to a method of class (1), the Schwarz-Christoffel 
method of L.N. Trefethen [DrTrj and its MATLAB implementation, the SC Toolbox 
written by T. DriscoU [Drij . 

In the two test cases we have used, the performance of the SC Toolbox was 
superior to AFEM. On the other hand, the AFEM software applies also to compu- 
tation of moduli of polygonal ring domains as shown in |BSVi| . AFEM also has an 
advantage in the problems where the quadrilateral has a large number of vertices. 
This situation arises when approximating nonpolygonal quadrilaterals (e.g. Example 
13.71) . We will report our results also in |RV2] . 
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2. Preliminaries 



A Jordan domain D in C with marked (positively ordered) points 2:1, 2:2, ^3, ^4 
e dD is a quadrilateral and denoted hy {D; zi, Z2, z^, z^) . We use the canonical map 
of this quadrilateral onto a rectangle (-D'; 1 + ih,ih, 0,1), with the vertices corre- 
sponding, to define the modulus h of a. quadrilateral {D; zi, Z2, z^, z^) . The modulus 
of {D; Z2, Z'i, Zi, zi) isl/h. 

We mainly study the situation where the boundary of D consists of the polygonal 
line segments through 2:1,2:2,^3,^4 (always positively oriented). In this case, the 
modulus is denoted by QM(D; zi, Z2, -23, z^). If the boundary of D consists of straight 
lines connecting the given boundary points, we omit the domain D and denote the 
corresponding modulus simply by QM(2:i, 2:2, ^3, Z4). 

The following problem is known as the Dirichlet-Neumann problem. Let D be 
a region in the complex plane whose boundary dD consists of a finite number of 
regular Jordan curves, so that at every point, except possibly at finitely many points 
of the boundary, a normal is defined. Let ■?/' to be a real-valued continuous function 
defined on dD. Let dD = AU B where A, B both are unions of Jordan arcs. Find a 
function u satisfying the following conditions: 

(1) M is continuous and differentiable in D. 

(2) u{t) = ij{t), t e A. 

(3) If d/dn denotes differentiation in the direction of the exterior normal, then 



One can express the modulus of a quadrilateral {D; zi, Z2, z^,, z^) in terms of 
the solution of the Dirichlet-Neumann problem as follows. Let 7j,j = 1,2,3,4 be 
the arcs of dD between (2:1,^2), (2:2,^3), (2:3,24), (^4,^1), respectively. If u is the 
(unique) harmonic solution of the Dirichlet-Neumann problem with boundary values 



of u equal to on 72, equal to 1 on 74 and with Sw/Sn = on 71 U 73 , then by [Ah\\ 
p. 65/Thm 4.5]: 



We also have the following connection to the modulus curve family (see e.g. [AVV^ 
pp. 158-165]): QM^D; zi, Z2, Z3, Z4) = M(r), where F is the family of all curves 
joining 72 and 74 in D. 

Another approach is to use the Schwarz-Christojf el formula to approximate the 
conformal mapping / onto the canonical rectangle. This formula gives an expression 
for a conformal map from the upper half-plane onto the interior of a n-gon. Its vertices 
are denoted wi, . . . , Wn, and aivr, . . . , a^vr are the corresponding interior angles. The 
preimages of the vertices (prevertices) are denoted hj zi < Z2 < . . . < Zn- The 



d_ 

dn 



t e B. 



(2.1) 




EXPERIMENTS WITH MODULI OF QUADRILATERALS II 



3 



Schwarz-Christoffel formula for the map / is 

where c is a (complex) constant. The main difficulty in applying this formula is that 
the prevertices Zj cannot, in general, be solved for analytically. By using a Mobius 
transformation, one may choose three of the prevertices arbitrarily. The remaining 
n— 3 prevertices are then obtained by solving a system of nonlinear equations. Several 
methods for solving this problem are discussed in |Bis] . |DrTr] . and |DrVa] . 

The MATLAB toolbox by T. DriscoU |Dri] contains a collection of algorithms 
for constructing Schwarz-Christoffel maps and computing the moduli of polygonal 
quadrilaterals. The toolbox also gives an estimate for the accuracy of the numerical 
approximation of the modulus. 



3. Experiments 

The solutions of the Dirichlet and the Dirichlet- Neumann problems can be ap- 
proximated by the method of finite elements, see [Hen[ pp. 305-314], [Pap]. Hence, 
this method can also be used to approximate the modulus of quadrilaterals and rings. 
The Dirichlet-Neumann problem can be numerically solved with AFEM (Adaptive 
FEM) numerical PDE analysis package by Klas Samuelsson. This software applies, 
e.g., to compute the modulus (capacity) of a bounded ring whose boundary com- 
ponents are broken lines. Examples and applications for this software are given in 
[BSV] . 

In |HVVj a theoretical formula for computing QM(A, S, 0, 1) was given with its 
implementation with Mathematica. This led to a study of the modulus of quadri- 
lateral in |DVj . In the course of the work on |DVj . the variation of the modulus 
was studied when one of the vertices varies and others are kept fixed, and several 
conjectures were formulated. For these purposes, neither the theoretical algorithm 
in [HVVj nor the implemented Mathematica program based on it was no longer ad- 
equate and we started to look for a robust program to compute QM(y4, B, 0, 1) . It 
seems that the AFEM software of Samuelsson is very efficient for this purpose. As in 
[RVlj we use the AFEM software of Samuelsson for computations involving moduli 
of polygonal quadrilaterals. 

3.1. Example. Let /(x, ?/) = QM(x+iy, i, 0, 1) - l/QM(?/ + ix, z, 0, 1). Then by [Hinl 
p. 433] we see that f{x,y) = 0. Therefore we may use this function as a measure of 
the accuracy of AFEM software and SC Toolbox. 
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Figure -659: Function logio(|/(x,t/)| + 10-^°) for x G (0,3], (0, 3] with AFEM 
(left) and SC Toolbox (right). 

3.2. Example. We study the function 

g{t,h) = QM{l + he'\he'\0,l). 
An analytic expression for this function has been given in |AQ VV , 2.3]: 

(3.3) git,h) = X'{n/^)/%{n/^), 

where 

7ih 



(3.4) 



, for < a < 1/2, 



2 sin(7ra) 

and the decreasing homeomorphism /i^: (0, 1) (0, oo) is defined by 

71 F(a, l-a;l;l-r2) 



(3.5) 
Here 



^"^""^ ~ 2sin(7ra) F{a,l - a;l;r^) 

oo 

F(a, b; c; z) = 2-^1(0, b; c;z) = '^ 



n=0 



a, n)(b, n) z"' 
(c, n) nl 



is the Gaussian hypergeometric function, 

(a, n) = a{a + l)(a + 2) . . . (a + n — 1), (a, 0) = 1 for a 7^ 0, 
is the shifted factorial function, and the elliptic integrals 'X{r),%'{r) are defined by 



3<:(r) = -F(l/2,l/2;l;r2), X'{r) = X{r'), andr' = Vl^. 
2 
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Figure -658: Function /ia(r). 

The function g{t,h) is the modulus of the parallelogram with opposite sides 1 
and h, respectively, and we see that there are three cases h G (0, 1), h = 1 and h > 1. 
In the first case the function is monotone increasing with respect to t G (0,7r/2), in 
the second case the function g{t, 1) = 1 is constant and in the third case decreasing. 




Figure -657: Function g(t, h) for t G (0, 7r/2) and h G [1/2, 2] (left), and the error 
estimate logio(|5'exact(i, h) - gnnmerit, h) \ + 10~^°) for the function g{t, h) (right). 
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4.5 




Figure -656: Function g{t, 1.5) for t G (0, tt/2). 



3.6. Example. The modulus QM(1 + i\a — 1|,2|6|,0,1) has an analytic expression 
if |a — 1| = h = |6| + 1 . Bowman [Bowl pp. 103-104] (see also |Burl p. 197]) gives a 
formula for the conformal modulus of the quadrilateral with vertices 1 -|- hi, {h — 
0, and 1 when h > 1 as M{h) = DC(r)/DC(r') where 

r=^^j, ii = /^v42^j' ^^ = ^v4Yj' c = 2h-l. 

Therefore, the quadrilateral can be conformally mapped onto the rectangle l+iM{h), 
iM{h), 0, 1, with the vertices corresponding to each other. It is clear that h — 1 < 
M(h) < h . The formula 

M{h) =h + c + 0{e~^^), c = -l/2-log2/7r ^ -0.720636, 

is given in [PSj . As far as we know there is neither an explicit nor asymptotic formula 
for the case when the angle 7r/4 of the trapezoid is equal to a G (0, 7r/2) . We compute 
the modulus QM(i/i, i{h — l), 0, 1) by using the square frame capacity formula, AFEM 
and Schwarz-Christoffel Toolbox. 
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h 


AFEM 


SC 


Accuracy 
/SC 


Bowman 


Error 
/AFEM 


Error/SC 


1.1 


0.3403159 


0.3403135 


1.787 


10- 


8 


0.3403135 


2.41655 • 


10- 


6 


1.57002 


10- 


9 


1.2 


0.4614938 


0.4614926 


1.734 


10- 


8 


0.4614926 


1.20727- 


10- 


6 


2.98441 


10- 


-9 


1.3 


0.5704380 


0.5704374 


5.310 


10- 


8 


0.5704374 


5.83493 • 


10- 


7 


4.59896 


10- 


-9 


1.4 


0.6747519 


0.6747518 


1.046 


10- 


7 


0.6747518 


8.83554 • 


10- 


8 


6.24458 


10- 


-9 


1.5 


0.7769433 


0.7769434 


2.408 


10- 


8 


0.7769434 


1.10607- 


10- 


7 


3.39673 


10- 


-9 


1.6 


0.8780836 


0.8780838 


1.920 


10- 


9 


0.8780838 


1.53305- 


10- 


7 


8.10543 


10- 


-10 


1.7 


0.9786840 


0.9786842 


5.439 


10- 


10 


0.9786842 


2.41392 - 


10- 


7 


2.02109 


10- 


-10 


1.8 


1.0790020 


1.0790024 


2.102 


10- 


10 


1.0790024 


4.03325 - 


10- 


7 


4.94438 


10- 


-11 


1.9 


1.1791710 


1.1791715 


6.225 


10- 


11 


1.1791715 


5.22481 - 


10- 


7 


1.20739 


10- 


11 


2.0 


1.2792610 


1.2792616 


1.536 


10- 


11 


1.2792616 


5.71171 - 


10- 


7 


2.97451 


10- 


12 



Table 1 : Error estimates for AFEM and SC Toolbox with the square frame capacity 
formula. The accuracy of the estimate given by SC Toolbox is also consistent with 
the experiment. 

3.7. Example. Let Q be the quadrilateral whose sides are defined by two circular 
arcs in the upper and lower half plane, perpendicular to the unit circle at the points 
giS^ gi(7r-fl)^ gi(6»-7r)^ g-i6i ^gjj ^Yiq two circular arcs through r, -i, —i and — r, 

i, — see Figure 6. If a, 6, c, d are the points of intersection of these four circular 
arcs in the II, III, IV and I quadrants respectively, then Q = {Q; a, b, c, d) defines a 
quadrilateral in the unit disk with 

QM{Q;a,b,c,d) = (vr - 2/5)/p, 

where 

1 + u 2r 

p = 2 log , (3 = arccot -, and u = tan(0/2). 

1 — u 1 — 



,i(,-e) 



/ E 











4-- P- 






-r; 1 ; r 












^ 

-i 





Figure -655: The circular quadrilateral Q. 
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e 


AFEM 


Exact 


Error 


0.10 


7.592357 


7.597433 


5.076 ■ 


10- 


'■i 


0.15 


5.056044 


5.054357 


1.687- 


10" 


-3 


0.20 


3.784480 


3.779611 


4.869 ■ 


10" 


-3 


0.25 


3.010221 


3.012175 


1.954- 


10- 


-3 


0.30 


2.497983 


2.498368 


3.849 ■ 


10" 


-4 


0.35 


2.130426 


2.129465 


9.616- 


10" 


-4 


1.00 


0.620785 


0.620631 


1.531 - 


10- 


-4 


1.05 


0.575335 


0.575402 


6.645 - 


10" 


-5 


1.10 


0.533529 


0.533010 


5.185- 


10" 


-4 


1.15 


0.493114 


0.492934 


1.805- 


10" 


-4 


1.20 


0.454678 


0.454689 


1.114- 


10" 


-5 



Table 2: The modulus of the quadrilateral Q for r = 0.4. 
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